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ABSTRACT 

Spatial and temporal backscauer intensity (<j°) variations from ice growing 
on shallow lakes during winter 1991-92 near Barrow, NW Alaska, have 
been quantified for the first time using ERS-t C-band SAR data acquired at 
the Alaska SAR Facility. A field and laboratory validation program, 
including measurements of the thickness and structure-stratigraphy of the 
ice, indicates that <?° values are strongly dependent on whether the ice 
freezes to the lake bottom, or remains afloat Backscauer intensity decreases 
significantly when the ice grounds on the bottom. Strong backscauer from 
floating ice is attributed to a specular ice-water interface and vertically 
oriented tubular bubbles. During the spring thaw, backscauer undergoes a 
reversal; o° values from ice that was grounded increase, while o° values 
from ice that was afloat decrease. This phenomenon has not previously been 
reported. 
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1. INTRODUCTION 

The X-band and L-band radar signatures of ice-covered lakes in late winter 
on the Alaskan North Slope were first investigated over a decade ago using 
airborne real aperture radar. The radar images showed interesting spatial 
variations in radar returns that were attributed to the effects of bubbles and to 
whether the ice was floating or frozen to the lake bed (Sellmann cl al., 1975; 
Etachi el al., 1976: Weeks et al.. 1977, 1978, 1981). A multi-temporal series 
of X-band radar images also clearly demonstrated how the radar returns 
changed during the course of the ice growth season as the area of ice frozen 
to the lake bed increased (Mellor, 1982). 

The previous investigations of the radar returns from the tundra 
lakes relied exclusively on uncalibrated photographic products; thus, it was 
not possible to quantify the spatial and temporal variations in backscauer 
intensity, i.e., o° values. In January 1992 radiomctrically calibrated, digital 
ERS-1 SAR data were released from Alaska SAR Facility to the scientific 
community. SAR data acquired at ASF during winter 1991-92 included a 
time-series of scenes of the Barrow region (Fig. 1) in north-west Alaska 
where there are hundreds of shallow tundra lakes. Spatial and temporal 
variations in backscauer intensity from ice growing on those lakes have been 
derived from ERS-1 digital SAR scenes. 

This paper presents the results of the field, ice characterization and 
SAR data analysis program. The ice characteristics and SAR data have also 
been used to formulate and test the first numerical model of backscauer from 
lake ice, which is described elsewhere in this conference proceedings 
(Wakabayashi et al., 1993). 

2. METHODS 

Field work was undertaken near Barrow in late April 1992. The purpose of 
the fieldwork was to locale the floating and grounded portions of selected 
lakes and to obtain ice cores for analysis of ice crystalline structure and air 
inclusion variations in the floating and grounded ice. The depth and state of 
the snowcover was noted at each sampling location. A total of 11 sampling 
sites were visited on seven lakes (Fig. 1). The 1992 field data are 
supplemented by observations made at a number of lakes in spring 1990 
(Fig. 1). 

For each lake that was sampled in late April the backscauer 
intensity (<T°) was determined from seventeen low-resolution (100m pixel 
size, 240m resolution) scenes obtained during winter 1991-92, beginning in 
late September 1991 when the lake ice cover was first established. A time- 
series of digital number (DN) profiles across Sukok Lake (Fig. 1) was also 
compiled to illustrate spatial changes in backscauer during the ice growth 
season. DN values are proportional to o° values. 



Fig. 1. Location map of sampling sites on lakes near Bartow, Alaska, in 
April 1992 (dots) and May 1990 (crosses). 

3. LAKE ICE STRUCTURAL-STRATIGRAPHIC CHARACTERISTICS 
3.1 Ice and snow characteristics 

Variations in snow depth, ice thickness and whether the ice was floating or 
grounded on the lake bottom at the lime of sampling are summarized in 
Table I. The ice was floating at six sampling sites and was frozen to the lake 
bottom at five others. Where water was present below the ice it occurred as 
only a very thin layer, sometimes only a few millimetres in thickness. The 
snowcover was dry and packed hard in most cases. Occasionally, the wind 
had blown the snow completely away, exposing patches of bare ice. 


Table 1. Summary of snow and ice characteristics at each sampling site 


Sampling 

Sampling 

Snow depth 

Ice thickness 

Floating or 

site 

date 

(mm) 

(m) 

grounded 

1 

19 April 

30 

1.945 

Floating 

2 

19 April 

60-80 

1.59 

Floating 

3 

19 April 

30 

0.53 

Grounded 

4 

19 April 

26 

0.24 

Grounded 

5 

20 April 

0-50 

2.1 

Floating 

6 

20 April 

80 

1.18 

Grounded 

7 

20 April 

20 

0.52 

Grounded 

8 

20 April 

50 

2.005 

Floating 

9 

28 April 

25 

2.16 

Floating 

10 

29 April 

60 

2.04 

Floating 

- 11 

29 April 

30 

0.77 

Grounded 


3.2 Crystalline Structure and Stratigraphy 

The ice cores contained three particular types of ice: 

1) granular ice containing roughly spherical bubbles, probably snow-ice; 

2) inclusion-free ice, some with the columnar texture indicating congelation 

crystals with horizontal c-axes, and some with the massive, single 
crystal structure typical of congelation ice with vertically-oriented 
c-axes; 

3) columnar and massive single congelation ice crystals containing tubular 

(cylindrical) bubbles oriented roughly normal to the ice surface, i.e. 
parallel to the ice growth direction. 
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The granular ice occurred at the top of cores from six locations 
(sites 2, 3, 5. 6, 7 and 11, Fig. 2). Clear ice occurred at all sites, either as the 
uppe r most layer or directly below the bubbly, granular ice (Fig. 2). Ice with 
tubular bubbles occurred as the lowermost layer at all locations (Fig. 2). 

The thickness of the granular layers varied between 0.03m and 
0.22m. The inclusion-free layer thicknesses ranged from 0.14m to 1.18m. 
The depth below the surface at which tubular bubbles appeared varied from 
0.26m to 1.18m. Tubular bubbles appeared earliest at the shallowest lake 
sites where the ice was frozen to the bottom when the samples were obtained 
(Fig. 2). The thickness of the ice layers with tubular bubbles ranged from 
0.27m to 132m. 

33 Bubble characteristics 

In the surface granular layers, the mean bubble diameter ranged from 0.7mm 
(site 2) to 4.6mm (site 3), and the number of bubbles varied between 
10.8 x lO^-m^ to 2,471 x Kp-m^. The density of tubular bubbles at the base 
of the ice varied between 0.4 x tO^-m^ to 41 1 x l(P-n|2. The radius of the 
tubular bubbles ranged from 0.13mm to 1.2mm. The largest diameter 
bubbles were found only at site 2 (Fig. 3); otherwise, the majority of tubular 
bubbles had a small radius averaging 0.4mm. The mean length of tubular 
bubbles ranged bom 140mm at Site 1 1 to 91mm at Site 2. The shortest and 
longest bubbles had lengths of 7.0mm and 178.0mm respectively. At any 
given sampling site, the density and mean diameter of the tubular bubbles 
varied with depth (Fig. 3). Consequently, the ice with tubular bubbles 
comprised a number of sub-layers (Fig. 2). 

4. VARIATIONS IN BACKSCATTER INTENSITY 

4.1 Temporal changes in backscatler intensity (o°) 

For each lake that was sampled in late April, changes in backscatter intensity 
as a function of time were determined for the central portion of the lake. The 
results for selected lakes are shown in Fig. 4. 

In late September when ice began to grow on the lakes, o° values 
were in the range -16 to -22 dB. By mid- to late October, values at all sites 
had increased to -12 to -15 dB. From lace October onward, o° values at site 
7 decreased further before levelling off in late December at -18 dB, where 
they remained throughout the winter. At site 11, maximum <j° values of 
-12 to -13 dB were reached in early December, then decreased during 
December and remained constant at -17 to -18 dB from early January to late 
April. 


on sitei Suez arras site io 




Fig. 2. Stratigraphy diagrams for ice cores illustrating the amount and 
distribution of ice types present, and those locations where the icc was 
frozen to the lake bottom. The horizontal lines in the tubular bubble layers 
represent sub-layers within that ice. 
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Fig. 3. Profiles of bubble density and bubble diameter at selected sites. BE is 
the mean length of the tubular bubbles. Sites 8 and 9 were on the same lake 
(see Fig. 1). 


Backscatter intensity at site 6 reached a maximum of -5 dB in early 
January, but by early February had decreased to -17 dB where it remained 
quite constant through late April. At Sites 1 and 2, o° values had increased 
to -5 to -6 dB by late January and remained quite constant at those values for 
the rest of the winter. 

From late April to early June, backscatter intensity went through a 
reversal. In late May, backscatter intensity at sites 6, 7 and 1 1 had risen to 
—12/— 14 dB, while at sites 1 and 2 it had fallen to —12/— 14 dB. By early June, 
backscatter intensity at sites 6, 7 and 11 had reached -7/-10 dB, while at 
sites 1 and 2 it had decreased further to -17/-19 dB. 

4.2 Spatial and temporal backscatter change at Sukok Lake 

The time sequence of spatial changes in digital number (DN) values for 
Sukok Lake shown in Fig. 3 provides a semi-quantitative illustration of a 
number of the items noted in the previous section. In late September the DN 
values are very low (minimum 40) across most of the lake, but as the autumn 
progresses they rise and by 1 December have increased to 140-150 across 
most of the lake (Fig. 5a). 

Between the beginning and the end of December there is a 
significant increase in DN values, from 140-150 to 180 respectively, but at 
the same time DN values, at the lake margins drop to 100-140 (Fig. 5b). 
From late December to late January there is a steady increase in DN values 
to as much as 200-210 across the lake, with a strong discontinuity at the 
margins where the DN values remain at 100-140 (Fig. 5b). 

From the end of January to late April there are significant changes 
in the DN values (Figs 5c and 5d). The marginal area of low DN values 
increases as the DN values fall to as low as 85, while the central area of high 
DN values shrinks as the DN values decrease to 175. 

4.3 Ice growth, structure and backscatter variations 

At each site that was sampled in 1990 and 1992, ice that was frozen to the 
lake bed corresponded to areas of low backscatter in the SAR scenes, and ice 
that was floating corresponded to areas of high backscatter. These 
observations add to the evidence that the low backscatter is due to the low 
dielectric contrast of the ice-soil interface, resulting in strong attenuation of 
the radar signal and minimal returns to die radar antenna (Scilmann ct al., 
1975; Elachi et al., 1976; Weeks et al., 1977, 1978, 1981; Mellor, 1982). 
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First day of mean daily temperatures above 0°C 
Fig. 4. Backscaaer coefficient (o°) changes as a function of time for a selection of lakes. 
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Fig. 5. Profiles of spatial changes in digital number (DN) values along the same transect across Sukok Lake 
frwn late September 1991 to 20 April 1992. The Julian dates (JD) correspond to the following calender 
dates: 272:29 September 1991; 290:17 October 1991; 299:26 October 1991; 335: 1 December 1991; 364, 30 
December 1991; 14:14 January 1992; 29:29 January 1992; 56:25 February 1992; 71:11 March 1992; 89:29 
March 1992; 1 1 1:20 April 1992. 
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Backscatter from the ice was also low at the beginning of the 
growth season (Figs 4 and 5).This can be attributed to specular reflection 
away from the rad ar antenna off the smooth horizontal surfaces (upper and 
lower) of the thin ice. As the ice thickened, backscatter values rose steadily, 
except for those lakes that were subsequently found to be frozen to the 
bottom. Backscatter from those lakes fell at different times during the winter 
(Fig. 4) when, it can be safely assumed, the ice first froze to the bottom. The 
changing backscatter profiles across Sukok Lake (Fig. 5) can be interpreted 
as a record of the progressive freezing of the ice to the lake bottom. 

At all sites where the ice was frozen to the lake bed, tubular bubbles 
were observed in the ice. This indicates that the tubular bubbles serve largely 
as forward scatterers as the radar signal passes through to the base of the ice. 
Since there is a minimal return to the radar once the ice is frozen to the 
bottom, the strong backscatter that is observed in the ice prior to grounding 
and in the ice that remains afloat all winter must require the presence of both 
an ice- water interface and the tubular bubbles. 

An ice-water interface has a strong dielectric contrast that will give 
rise to specular reflection away from the radar (Weeks et al„ 1978). 
However, at the time the lakes were sampled, there was strong backscatter at 
those sites where the ice was still floating. The only explanation for the 
strong backscatter is the presence of the tubular bubbles, which would cause 
volume scattering of the radar signal once it had been reflected off the ice- 
water interface (Weeks et a)., 1978). 

The significant role of the tubular bubbles is illustrated by the 
backscatter differences between the ice that remained afloat and that which 
was frozen to the lake bed at the end of April. With o° values of -5 to -6 dB, 
backscatter from the floating ice is significantly greater than the backscatter 
from the grounded ice (-15 to -18 dB). A similar large difference exists 
between the backscatter from ice that remained afloat in April and 
backscatter from the same ice at the beginning of the growth season (Sites 1 
and 2, Fig. 4). It is also worth noting that the backscatter from the grounded 
ice is similar to the backscatter from the new ice that began to grow in 
September (Fig. 4). This similarity in backscatter provides an indication of 
the effectiveness of the new ice as a specular reflector prior to tubular bubble 
development 

• Earlier studies of the SAR signatures of ice on shallow Arctic lakes 
have emphasized the importance of the tubular bubbles (Sellmann et al., 
1975; Elachi et al.. 1976; Weeks et al.. 1977, 1978, 1981; Mcllor, 1982). No 
previous mention has been made of the possible impact on backscatter 
variability of the bubbly snow-ice layer at the surface (Fig. 2). Some 
indication of the relative importance of the bubbly surface layer can be 
gained from a comparison of backscatter at Sites 1 and 2. Both sites were on 
the same lake (Fig. 1), but only Site 2 had a bubbly surface layer. Despite the 
physical difference between each site, there is a negligible difference in 
backscatter at each site. Thus, it appears that the effects of the tubular 
bubbles far exceed any effects due to the presence of bubbly snow-ice 
layers. 

Unusual changes in backscatter intensity have been recorded during 
the spring thaw period in May and June; i.e., backscatter from those lakes 
that were grounded for most of the winter increased, while backscatter from 
those lakes that had remained afloat all winter decreased. This backscatter 
reversal has not previously been observed and as yet we have no explanation 
for the phenomenon, other than that it probably relates to changes in the state 
of the ice. 

5. CONCLUSION 

A winter 1991-92 time-series of radiometrically calibrated ERS-1, C-band 
SAR data of the Barrow, Alaska region was acquired at the Alaska SAR 
Facility. Spatial and temporal backscatter intensity variations from ice 
growing on lakes near Barrow have been derived for the first time. The SAR 
data show strong backscatter variability related to ice growth, development 
of tubular bubbles and freezing to the bottom. There is strong attenuation of 
the radar signal once a lake freezes to the bottom, regardless of the presence 
of tubular bubbles in the ice. The tubular bubbles are effective scatterers 
providing a strong return signal to the radar only as long as the ice is floating 
and a specular ice-water interface is present Time-series of backscatter 
variations from these shallow tundra lakes arc, therefore, a record of, 1) the 
development of tubular bubbles in the ice and, by association, gas saturation 
of the underlying water, and 2) the freezing of ice to the bottoms of the lakes 
and, therefore, lake bathymetry. The availability of a radiomctrically- 
calibrated SAR time series has also allowed us to report for the first time an 
unusual backscatter reversal during the thaw period. The reversal is probably 
due to changes in the nature of the lake ice cover which remain to be 
determined. 
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